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^SJ . The effect of fluctuations about the tt-Aux mean field state for the undoped high- 

er . temperature superconductors is investigated. It is shown that fluctuations of the mean fields 

D ■ 
/^ , lead to a self-energy correction that doubles the band width of the fermion dispersion in the 

m 






lowest order. The dynamical mass generation is associated with the self-energy effect due to 
the interaction mediated by the Lagrange multiplier field, which is introduced to impose the 
constraint on the fermions. A self-consistent picture about the mass generation and the prop- 
agation of the Lagrange multiplier field without damping is proposed. The antiferromagnetic 
^ , long-range ordering is described without introducing an additional repulsive interaction. The 

t^ , theory suggests a natural framework to study spin disordered systems in which fermionic 

excitations are low-lying excitations. 
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1. Introduction 



(N 
> 

in 

en . In the phase diagram of the high-temperature superconductors, apparently the most es- 

cn ' tablished phases are the d-wave superconducting phase and the Neel ordering phase in the 

^r I undoped compound. One of the key questions about high-temperature superconductivity is 

00 , how to connect these two phases. If one starts from the Neel ordering state toward the d-wave 

O ■ 

superconducting state, the first step would be to consider the single hole doped system. 

Experimentally such the system has been studied by angle-resolved photoemission spec- 
H ' troscopy (ARPES) in the undoped compounds^"^ where a photo-hole is introduced in the 

system and the excitation spectrum is associated with the properties of the single hole doped 
system. The experimentally obtained excitation spectrum is in qualitatively good agreement 
with the quasiparticle excitation spectrum of the vr-flux mean field state^'^ with a mass term.^ 
Although at the mean field level the quasiparticles are gapless, an excitation energy gap opens 
up by adding an on-site repulsive interaction. The gap arises from the staggered magnetiza- 
tion, and by taking it as a variational parameter, better variational energy is obtained.^ An 
effective field theory approach also suggests the presence of the mass term.'^ In addition, the 
dispersions along (0, 0)-(7r, vr) line and (0, 7r)-(7r,0) line in the Brillouin zone are isotropic as 
observed in the experiments.^ (The lattice constant is taken as the unit of length throughout 
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the paper.) However, the band width of the quasiparticles in the vr-flux mean field state is 
smaller than the experiment by a factor of 2 — 3. Furthermore, there is no reliable estimation 
of the mass value except for variational wave function approaches. However, in the variational 
wave function approaches it is necessary to include an additional repulsive interaction. 

On the other hand, self-consistent Born approximation analysis of the t-J model*^ sug- 
gests that the band dispersion is significantly renormalized from t ~ 0.4eV to J ~ O.lSeV, 
which is consistent with the experiments. However, there is some discrepancy in the excita- 
tion spectrum along the line in the momentum space from (vr, 0) to (0, vr). This discrepancy is 
removed by including the next-nearest and the third nearest neighbor hopping processes.^ In 
this approach, the t-J model is analyzed in terms of the slave-fermion mean field theory. In the 
slave-fermion mean field theory, the spins are described by the Schwinger bosons. ^"^ Therefore, 
it is straightforward to describe the Neel ordering state as Bose-Einstein condensate of those 
bosons. However, to describe the d-wave superconducting state the slave-fermion formalism is 
not convenient. For the description of the d-wave superconducting state, slave-boson formal- 
ism is used in the literature. ^^ In order to avoid taking a different formalism, here I focus on 
the TT-fiux state. 

Before going into discussions about the vr-flux state, let us discuss advantages and disad- 
vantages of the theory. In order to consider the single hole starting from the spin 1/2 anti- 
ferromagnetic Heisenberg model, mainly there are two approaches: One is to represent spins 
in terms of boson fields, like in the non-linear a model^^ or equivalently in the CP^ model^^ 
derived from the Schwinger boson mean field theory.^'' In this approach, we need to think of a 
soliton-like excitation to describe a doped hole, which is fermion, in terms of boson fields. The 
situation is very similar to the Skyrme model^^ in which fermions (protons and neutrons) are 
described as solitons in the non-linear a model. (Goldstone modes, which correspond to spin 
wave excitations in the Heisenberg antiferromagnet, are pion fields.) It is argued that in ref. 
15 that a doped hole can be described by a skyrmion-like spin texture. ^^"^^ The band width 
and the mass, which is the excitation energy of the spin texture, are in good agreement with 
experimentally estimated values. (In the experiments, the mass value can be estimated by 
approximating the dispersion around (ib7r/2, ±7r/2) by a conventional non-relativistic kinetic 
energy form.) From the Ginzburg-Landau description of the Heisenberg antiferromagnet, it 
is natural to expect the appearance of such a spin texture in the continuum: In the bosonic 
theory the Neel ordering state is described by a Bose-Einstein condensate. It turns out that 
local suppression of the order parameter, which is introduced by the formation of the Zhang- 
Rice singlet, ^^ leads to a vortex-like state. This solution is easily found from the analysis of 
two-dimensional Ginzburg-Landau theory. ^^ From the CP^ formalism, ^^ the vortex like state 
turns out to be a skyrmion like spin texture. However, stability of such a state on the lat- 
tice is not evident. This is because that the bosonic theory describes mainly the low- lying 
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excitations, that is, the spin wave excitations. Therefore, it is not easy to study microscopic 
stabiUty of the spin texture within a bosonic theory. 

By contrast, the other approach is to describe spins in terms of fermions: 

Sja = ^fj'^afj-, (1) 

where /• = (/Jt)/J|)) ^^a (a = x,y,z) are the Pauh matrices. Because the spins are 1/2, 
there is a constraint on the fermions. 

It is expected that the quantum nature of the spin 1/2 is well described by this theory. How- 
ever, there is disadvantage in the description of the antiferromagnetic long-range ordering. 
Within the mean field theory, the staggered magnetization vanishes. Therefore, we need to 
go beyond the mean field theory to describe Neel ordering. In the continuum limit, the Neel 
ordering is discussed^ in the context of the dynamical mass generation in quantum electrody- 
namics in three spatial and time dimensions. ^^"^^ Although there are discussions about the 
condition of the mass generation in the sense of the 1/N expansion, it is hard to estimate 
the value of the dynamically generated mass. Another way to describe the mass generation is 
to include the short-range Coulomb repulsion.^ However, it is not clear whether adding the 
short-range Coulomb repulsion term to the antiferromagnetic Heisenberg model is necessary. 
The antiferromagnetic Heisenberg model is derived from the Hubbard model taking the limit 
in which the on-site Coulomb repulsion, U, is much larger than the hopping matrix element, t. 
The condition of U/t ^ 1, is replaced by the constraint. So, it is unclear whether an additional 
repulsive interaction is necessary. 

In this paper, we study the fluctuation effects in the vr-flux state. The constraint is included 
in the action of the system using a standard Lagrange multiplier field. The effect of the 
constraint is studied by analyzing the self-energy effect associated with fluctuations of the 
Lagrange multiplier fields. It is argued that by including fluctuations the resulting quasiparticle 
band width and the staggered magnetization are in good agreement with the experiment. 
The result suggests that the vr-flux state is a promising candidate for the description of the 
undoped high-temperature superconductors. In particular, the fact that the model described 
by the fermion flelds includes the Neel ordering state implies that the theory provides us a 
natural framework to study doping effect and the spin disordered phase in the presence of the 
doped holes. Because the quantum nature of the holes are well described by such the theory. 

The organization of this paper is as follows. In Sec. 2 we study the fluctuations of the 
vr-flux state mean flelds. It is shown that the lowest order self-energy correction doubles the 
amplitude of the mean fields. In Sec. 3 the coherent state path-integral action is introduced as 
a systematic approach. The effect of the Lagrange multiplier field and the vr-flux state mean 
fields is studied in Sec. 4, and the dynamical mass generation is discussed. In Sec. 5 the spin 
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wave excitation, which is associated with phase fluctuations of the vr-flux state mean fields, is 
discussed. Finally, section 6 is devoted to summary and discussion. 

2. The TT-Flux State and Perturbative Analysis of Fluctuations 

In order to investigate the effect of fluctuations, let us start with the Hamiltonian of the vr- 
flux state. Here we derive the mean field Hamiltonian and the term that describes fluctuations 
about the mean field state. The spin 1/2 antiferromagnetic Heisenberg model on the square 
lattice is given by 

H = JY,S^■SJ, (3) 

where J is the superexchange interaction between the spins. The summation is taken over the 
nearest neabor sites. As described in Introduction, the spins are represented by the fermions 
defined by eq.(l) with the constraint (2). In terms of these fermions, the Hamiltonian reads, 

^ = -^-^E E flh^flhp + E ^^ ( E fU=^ - 1 ) ' (4) 

up to a constant term. Here the constraint (2) is included by using a Lagrange multiplier. 
Now we introduce the following mean fields, 

(J 

where S = ±x, ±y. The vr-flux state is obtained by assuming 

Xj,j+x = Xj,j-x = X, (6) 

Xj,j+y = Xj,j-y = iX, (7) 

and Xj =0. After the Fourier transform the mean field Hamiltonian is given by 

HmF = E E( fka flk,a ) _ * n i 

keRBZ a \ '^k ^ / 

x( ^:'-]+NJ\x\\ (8) 

\ Joka / 

where the summation with respect to k is taken over the reduced Brillouin zone, \kx — ky\ < vr 

and \kx + ky\ < vr, and 

f^k = xJ (cos kx + i cos ky) . (9) 

The fields feka and foka are defined as 

feio)k. = J^ E /.^e"^'■''^ (10) 

j€A{B) 
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where the two sublattices are labeled by A and B. The mean field equation is given by 

1 = - y ^^!!^iL±^^!%tanh^. (11) 

keRBZ 

At T = 0, we find x — 0.479. As discussed in Introduction, the fermion dispersion band width 
using this mean field value is half of the experimentally estimated value. 

In order to calculate the self-energy correction associated with fluctuations, we study the 
Green's function. Introducing two component spinor, 

Jka = y feka f oka ) ' \^'^> 

we define the Matsubara Green's function as follows, 

Gka{r) = -{fka{r)flm. (13) 

In the mean field state the Green's function is 



^fcT(^^")= L_,*r 1 1 ^ M , ni ' I' (14) 



^ k \ iWn+Ek iWn-Ek ) iuln+Ei^ ' iidn-Ef; 



„(0) /• \ I iuin+Ek ^ iu}„-Ek '^ ^ \iuj„+Ek iuj„-Ek i i /-. r-x 

^kli^^n)= \ ^_^^( ^ ^ ^ |,\|2 ^ „2 ^ ), (15) 



fc k yiuin+E,, ~ iuJn-Ek ) iuj„+Ek ' 

where Uk = l/\/2 and v^ = Kk/{\/2\Kk\)- 

Having defined the mean field Green's function, let us derive the term describing fluctua- 
tions. Fluctuations about the mean field, x^ arise from 

k,k',qeRBZ,q=/=0 
^Jo,k+q,aJokl3jf,kijiJek'+qa- (16) 

Note that g = contribution is excluded. Because the q = Q term is included in the mean field 
Hamiltonian. In order to make the notation simpler, we omit g 7^ in the following equations. 
The interaction vertices are represented in Fig. 1. For example, the diagram with a = TqCt^ 
and b = TeCr- in Fig.l corresponds to the following term, 

H(roa+,r,a-) = Jj Y^ {cOS q^ + COS Qy) 

k,k',qeRBZ 

^fo,k+q,^fokifek'ifek'+q't- (17) 

The importance of including the fluctuation effect is clearly demonstrated by the lowest 
order self-energy correction. The first order self-energy with respect to Hi^t for the spin-up 
fermions is 

Sfct (^t^n) = -;a7V S I]I]2Jq \Gk+qa{'i^n + i^n) [ 

qeRBZiUr. a \ 1 J \ J 

/ \ / \ ' ^ ^ 

q(iRBZiQ.r, o- \ / V 1 / 
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Fig. 1. The Fcynman diagram for the interaction term associated with fluctuations of Xj,j+s- Here 

(a, 6) = (To<T^,Te<T'f),{To<Tl,Te(Jl),{To(J^,Te<T+),{To(J+,Te<T^). 



Straightforward calculation leads to 

2 



Sfc| (iiVn) 



N 



E ^^ 



k'eRBZ 



k'-k 







2uk'Vk' \ 

2uk'vl, J 

—CJ cos kx — iCJ cos ky 



-CJcoskx + iC J cos ky 







(19) 



with C ~ 0.479. Thus, the band width of the quasiparticle dispersion, which is given by 



E^^ = xJy^cos^ kx + cos^ ky, is exactly doubled by this self-energy. Therefore, the quantita- 
tive difference of the quasiparticle band width compared to the experiment is much improved 
by including the fluctuation effect. 

3. Coherent State Path-Integral Formalism 

In principle, one can extend perturbative analysis in the previous section to higher order 
trerms. However, since there are eight vertices, drawing Feynman diagrams is complicated. 
Furthermore, the Hamiltonian formalism is not useful for the study of fluctuations of the 
Lagrange multiplier field. For this purpose, the coherent state path-integral formulation is 
convenient. In the coherent state path-integral formulation, the partition function is given by 

E = JDpDfDXexpl-S], (20) 

where the action is 



S 



13 



dT 



Ti/S 



(21) 



7 . fja (^r + Aj) fja - ^'J /^ fiafjafjpfil. 
3 (id) 

with (3 = l/ksT the inverse temperature. After the Fourier transform, the Lagrangian is given 

by 



ft 



ft 



ft 



ft 



^— Z^ [jeka^TJeka + Joka^'^Joku) N l^ Z^ Jok+qaJekaJek'pJok'+ql3^ 

k&RBZ ^ ' k,k',qeRBZS=±x,±y 

+ 1^ [ ^gJek+qaJeka- + \+QJek+qaJ'^ka + \jok+qaJok(7 " ^q+QJok+qaJoka 
k,q&RBZ ^ 



-i{k-k')-5 



(22) 
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where Q = (tt, vr). 

Now we introduce a Storatonovich-Hubbard transformation to rewrite the interaction term 
as follows, 

L = 2^ [feku^rfeka + foka^rfokuj 

keRBZ 

+ 2_^ 2_^ y\jek+qaJekcT + \+QJek+qaJeka + \jok+qaJoku — \+QJok+qaJoka 

qeRBZ keRBZ 



J_ 

2 



/ > 2^ ^9 [y ]\f Z^ ^ JekaJok+qa I 

qeRBZ S=±x,±y \ keRBZ / 

~2" ^ -^ '^'^ iViV ^ '^ fok+qafekaj 

qeRBZ S=±x,±y \ keRBZ / 

qeRBZ S=±x,±y 

At the saddle point, we see that 

"2 



keRBZ 

As a mean field state, we assume that 



4'^ = ViV ^ ^*'(/eL/o;^+.«)- (23) 



Xj"''-' = Y yX3<^g,0, Xj"^^ = \lYX4Sq,0, 

A, = 0, Ag+Q = 0. (24) 



The action is 



S = dT 

Jo 



t W ^^ '^fe 1 ( -^efc" 



/ ^ y •' efco J oka y i ft I \ f 

keRBZ \ l^k C'r / \ /ofca 



A' <V /, |2 I i2 I i2 I i2 

-^ — T- \\Xl\ +\X2\ +|X3| +\xA 



where 



(25) 



«fc = "2 (xie '^" + X2e *^^ + Xse*^" + X4e^^^ ) . (26) 



For the vr-flux state, Xj's are 

Xi = X3 = -iX2 = -iXA = X- (27) 

The saddle point equation, or the mean field equation, is given by eq. (fl). 
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4. Effect of Fluctuations 

Now we study the fluctuations of the Lagrange multipher field and the vr-flux state mean 
fields, X- Fluctuations about the saddle point are described by the following action, 



Sf = j^ dr 



2 2^ 2^ ^^^^^ "*" 2^ V ■'ek+qa- J ok+qa 

qeRBZ S=±x,±y k,q(^RBZ 

\ J_ 

S=±x,±y 



^ g-»(fc+f).5. 




o=±x^±y 

where Xq,& = e}'^'^''^Xq and Agg = Ag + A^+g and Aog = Ag — A^+q. We consider fluctuations 
within the Gaussian approximation. At this point, we include the staggered magnetization 
term: 

^0 ^^^S2 ^ ' \ ^ ^stCTz J \ fok J 

Within the mean field theory, Agt = 0. Nonzero value of Agt arises from the self-energy effect 
as shall be seen below. Integrating out the fermion fields leads to the following effective action. 



*/ - r2^[ ^-'J ^"1 J ttA /„a ttA 



where 






keRBZ V^^" "*" -^k+q + -Ek «^n - -E'k+q — E^ J ' 
keRBZ -^k-E'k+q \ i^n + -^k+q + -^k i^n — -E'k+q — E^, J ' 






^ keRBZ -^k-E'k+q V ^^n + -^k+q + -E'k i^n - -E'k+q " -E'k 

^ k&RBZ -^k-E'k+q \ ^^^n + -E'k+q + -E'k ^^n " -E'k+q " -E'k / 



Here £^k = y/x'^J'^ (cos^ /c^; + cos^ ky) + A^^, and g denotes (q, Ji7„) with q gRBZ and $!„ = 
2Trn/P being the bosonic Matsubara frequency. 

In order to make clear the physical meaning of Xq and XqS, we study them separately. 
First, let us investigate the effect of the former. As we will see below, exchange of Xq fields 
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Fig. 2. The inverse of the potential V ! ■, versus q/n. 



leads to a logarithmic confining potential between the fermions in the presence of the fermion 
mass term. The coupling term between Ag and the fermions is 



ST' 



E 

k,q 



Aro a, 



^ek+qcr/eko- 



eq 



^oq 



f, 




Jeko-^ek+qcr Joko-Jok 



-qa 



Ao,. 






8 



Aro a, 



veq 



og 



'ok+Qo-^oko- 

no"(g) -n^(g) 
-n^(g) n^(g) 

Note that using this action in deriving the effective interaction term is equivalent to the 
random phase approximation. This point is demonstrated using a simple model in Appendix. 
Integrating out Xq leads to the following interaction term, 




(36) 



Sx 



E 



where 



Peg ~r Poq Peq Poq 

Pe-q + Po-q \ 
Pe-q — Po-q I 



y(±) = I 

n^ {q) T nt (g) ■ 



Vc 



(+) 



Vc 





(-) 



(37) 



(38) 



0, the denominator of Vq behaves as 1/q^ for 



^(+) 



One can see that for the static case, i$7„ 

g ^ as seen from the expansion with respect to q. {1/Vq^^ is shown in Fig. 2.) Therefore, the 
potential term Vq is a logarithmic potential. ^^ Since such a potential diverges at long dis- 
tance, it is a kind of confining potential. The same instantaneous interaction term is obtained 
in the continuum theory of QED3 by taking the Coulomb gauge. 

Now we examine dynamical mass generation associated with this logarithmic confining 
potential. We consider the first order self-energy term. The equation to determine the self- 
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energy in a self-consistent manner is 

Sek {i(^n) = -^ ^ VJ^^Gek+a {i^n + i^n)- (39) 

q 

A similar equation holds for S^k (i^n)- Here GekiiuJn) = G^j^ (iu;„) — Sek(^'i^n)- To study the 
dynamical mass generation, we need to consider a non-uniform mass term in general. However, 
for such a general mass term we are unable to obtain analytic forms for the propagators. To 
solve the equation (39) self-consistently we need to solve the coupled equations numerically. 
Here we discuss the dynamical mass generation within the approximation of the uniform mass 
term, Sek(iw„) = —AgtO'z- Under this approximation, we find the equation for the mass, 



-T 



K 



(+) 



q 



^ q 2^(7r/2,7r/2)+q 



1. (40) 



In spite of this simple form, numerically solving this gap equation is not an easy task. Because 
the potential is singular, it is necessary to introduce an infrared cutoff and the result depends 
on it. Furthermore, the result also depends on the number of lattice points introduced for 
numerical estimations.^^ Instead of precisely determining the mass gap value, we evaluate it 
approximately. After some algebra, we find 

n^e(.)-n^o(.) = ^ E ^+(-'f + --'l)^ E ^' (41) 

where 

'ik,q - -^'Ja'J/3 y-Qk^QJ^^^^ ~W^Wf,, 

In the /c-summation in eq.(41), dominant contribution comes from k ~ (ib7r/2, zizir/2). Taking 
the approximate form around these points, we obtain 



N ^^ El \ 2 2J N ^^ Ei 



^k,q ^ -;<ic.qp 7j7-ij7-^k - u7-ur- • (42) 



1 ixJf 

2 A3 



Ke{q)-Ko{q)^7;'-^^- (43) 



Now the gap equation is 



2 

Si 

27r(xJ)2g 



^'^ ^'* -1. (44) 



Note that this equation is different from that for QED3 by the factor of (A^t/xJ)^ in which 
the potential term is not proportional to A^^. Apparently this gap equation is suffering from 
the infrared divergence. In order to properly deal with this infrared divergence, we need to 
take into account the vertex correction. In the QED3 theory there is controversy in the choice 
of the vertex correction. 3^' ^6 Here we determine it from a physical argument. The interaction 
term in the gap equation comes from the following term, 

H,^, = yvMsp,8p^,, (45) 
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where 5pq = pq — J), with p being the uniformely distributed background particle density. Note 
that 5pq -^ for q ^ 0. On the other hand, from the density-density correlation function, 
(dpqdp^q) ~ q'^. This suggests that the approximate form of the vertex correction is 



(A„/xJ)-°''<(" (for <,<v'a:;7Ij) 

''^^ 1 (for ,>vOT) • *'"' 



with a = 1. Dividing the integration with respect to q in eq.(44) with including 7^, into 
0<q< ^/lAst/xJ) and q > ^{^st/xJ), 

2,(A,„/x/)^-l„g(^^-^).l. (47) 

This equation can be solved numerically, and the solution is 

A,t ~ 0.616xJ. (48) 

Substituting % ~ 2 x 0.479, we find A^j ~ 0.59J. This value corresponds to the staggered 
magnetization oinist — 0.30, which is a reasonable value compared to the numerical simulation 
results and the experiments. 

Having discussed the dynamical mass generation associated with the Lagrange multiplier 
field fluctuations, let us move on to the mean field fluctuations, XqS- The fluctuation mode 
associated with Xq& is found from the pole of the following equation, 

/ l-n^(g) Ii}{q) \ , , 

det ° ^^^ ^ ^^' = 0. 49 

V n|(g) l-Iiliq) I 
The pole was found from numerical computations with mesh size of 64x64 and taking 5 = 0.10, 
where 5 is introdued in the analytic continuation from the Matsubara frequency, iQ.n to the 
real frequency, through iVtn = uj+i5. It was found that the amplitude fluctuation mode has the 
excitation energy gaps whose minimum is about 1.2 J. Along the line from (0, 0) to (vr/2, vr/2) 
in the reduced Brillouin zone, the band minimum is located at (0,0), and the maximum is 
at (7r/2,7r/2). The band width along this line is 1.3J. The excitation energy at other points 
are within this energy range. Although the amplitude fluctuation mode of XqS is high-energy 
mode and can be neglected, the phase fluctuation of XqS is related to the spin-wave excitation. 
Because there is the relation, 



(Sj • ^j+s) - IXj,i+<5| 



(l&RBZ 

The spin wave excitation is discussed in the next section. 



I ^ e-(^^+V2)^^^, 



(50) 



5. Spin- Wave Excitations 

In this section, we consider the spin wave excitation associated with the phase fluctuations 
of the mean field, XqS iii ths presence of the dynamically generated mass, or the staggered 
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magnetization. The calculation is simiar to the spin-density wave theory. '^^'^^ Here we calculate 
the spin wave excitation spectrum following ref.39. 

The repulsive interaction Vg that leads to the non-zero staggered magnetization is ap- 
proximated by a short-range repulsion, V, for simplicity. The interaction V is associated with 
a weak short-range repulsive interaction assumed by Hsu.^ The value of V is evaluated from 

1 V 1 _ 1 

keRBZ y JKkl + \Ast\ 
Using the values of x = 0.994 and A^^ = 0.616J, we obtain V = 2.2J. The action that 
describes fluctuations is given by 



S 



/3 



^^E fl f> 



t 

k+Q 



dr+xJ COS kx -ixJ COS ky- A stCTz \ ( fk 

ixJ COS ky - A stCz dr-x-J COS kx J \ fk+ 



I dT -^ Smk-k' ■ [fk^(yfk'j + 77 E '^™9 • <^ni- 



k,k' <? 

Integrating out fermions, and after some algebra, we obtain 



+ y ^Sf 



^eff 



q&RBZ 



+ - Yl ( ^^q ^1^l+Q ) 

qeRBZ \ 



-Cq bq+Q 



^^-q-Q 



(51) 



+ T7 Y. i ^"^q+Q ^"^9 



q&RBZ 



Cq bgj^Q 



^^-q-Q 

5m^ „ 



where 



with 



a„ = l--V 



q q 



,bg = l-\vK(^\cg = WK(-\ 



keRBZ \ '^^^ '^ / 



1 



1 



iUn + -E^k+q + -E'k i^n — -^k^ 



(52) 



(53) 



1 N Z^ 



'2A% 



N j^l^^ £^k+q^k \i^n + -E'k+q + -^k i^n - -^k+q - -E'k 



(54) 
(55) 

• (56) 
keRBZ 

The spin wave excitaion is associated with the pole of the transverse spin fluctuations. The 
pole is found from 

ba 






-iA 



St 



— ^—1 ( ^ 1 

Ek ^k+q / V ^^n + -^k+q + -E'k «^n — -E'k+q ~ -^k 



det 



-Cq bg^ 



0, 



(57) 
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Fig. 3. The mean field x dependence of the spin-wave velocity, Cs 



that is, 



i4^'^'™)0-H°'«) + (^"^t')' = o- 



(58) 



We perform the analytic continuuation of ii7„ 
respect to lo and q. Noting that 

2 Q 



uj + i6, and then expand each quantity with 



^K';:' {inn = 0) 



1 ^ 1 _ 1 



N ^^ El, 

keRBZ '^ 



v 



we find, after some algebra and a numerical computation, 

CO ~ 0.87 Jq. 



(59) 



(60) 



Thus, the spin-wave velocity is Cgw = 0.87J. This result is inagreement with the known 
established value of Cgw = 1-65J estimated from the quantum Monte Carlo simulations and a 
series expansion. '^'^"^^ However, the value of Cgw depends on x- Since the estimation of x based 
on a perturbative analysis, the value of x would change by further including higher order 
terms. Here we do not attempt to estimate x precisely in this way because it is hard to fix the 
value from that procedure. Instead, we compute x dependence of Csw by fixing Agt = 0.60J, 
which is the constraint from the experiments and the numerical simulations. In Fig. 3, x versus 
Csw is shown. From the known value of Cg^, x is evaluated instead. Cgw — 1-65 is obtained by 
setting X = 1-25. This value implies that the quasiparticle band width is 1.76J, which is a 
reasonable value compared to the experiment. 

In the limit of Agt/ixJ) ^ 1) the spin wave dispersion is easily obtained® as in the 
spin-density wave state. ^'^'^^ Using the approximation like 



1 



1 



\Kk\ 

2Af,' 
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Fig. 4. Numerically calculated spin wave dispersion along high-symmetry directions. The lattice size 
is 40 X 40 and S = 0.2. 



and noting 



1 V- * 

keRBZ 



Kk 



ixjy 



(cos Qx + COS Qy) , 



we find the following dispersion, 



UJ„ 



ixjy 



1 



1 - - {cosqx + cosqy) . 



(61) 



(62) 



This coincides with the result of the spin wave theory of the antiferromagnetic Heisenberg 
model. The same form of the dispersion is also obtained in the spin-density wave state in the 
strong coupling limit. ^^'^^ Although the assumption of Ast/{xJ) ^ 1 is not valid for excita- 
tions with to > Ast, eq.(62) is better for high-energy excitations than numerically obtained 
spin wave dispersion^ shown in Fig. 4 along high-symmetry directions in the Brillouin zone. 
The situation is similar to the spin density wave state as shown in Fig. 5. 

The full spin wave excitation spectrum is investigated by the neutron scattering experi- 
ments in the undoped compound. ^^ Most of features is in good agreement with the spin wave 
dispersion (62) with a suitable prefactor except around (7r,0). The dip around this point can 
be explained by including a ring exchange term to the Heisenberg model. ^^'^^ Effect of righ- 
exchange interaction on the vr-flux state is considered in ref.45 within the mean field theory. 
Although it would be interesting to investigate fluctuations about the mean field state includ- 
ing the ring-exchange interaction, consideration of such an effect is beyond the scope of this 
paper. 

6. Summary and Discussion 

In this paper, the effect of fluctuations about the vr-flux mean field state has been investi- 
gated. As for the fluctuations of the Lagrange multiplier field, it is shown that the fluctuations 
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Fig. 5. Numerically calculated spin wave dispersion along high-symmetry directions for the spin- 
density wave state at half- filling. The lattice size is 40 x 40 and 5 = 0.1. The on-site Coulomb 
repulsion \sU — ht. 

lead to a logarithmic potential term in the presence of the mass term. The mass value is eval- 
uated by solving the Schwinger-Dyson equation approximately, and a reasonable value is 
obtained. Because of the mass term, there is an energy gap in particle-hole excitations. As a 
result, the interaction mediated by the Lagrange multiplier field fluctuations is long-ranged. If 
the fermions were massless, the interaction would be short-ranged. (For a different approach 
to the constraint, see ref.46.) 

As for the fluctuations of the vr-flux state mean fields, the lowest order self-energy cor- 
rection doubles the mean field value. The resulting quasiparticle band withd is much closer 
to that evaluated in the ARPES experiment^ compared to the mean field result. From the 
pole of the correlation function, it is found that the amplitude fluctuation mode has a gap of 
1.2J. Meanwhile, the phase fluctuations are associated with the spin wave excitation, which 
is analyzed by approximating the potential term mediated by the Lagrange multiplier field 
fluctuations by a short-range repulsion. 

The dynamical mass generation is consistent with the confinement of fermionic excita- 
tions. That is, there are no low-lying fermionic excitations in the two-dimensional Heisenberg 
antiferromagnet. Due to the mass term of the fermions on the order of J, the propagation of 
the Lagrange multiplier field do not excite particle-hole pairs. Therefore, there is no damp- 
ing for the propagation. As a result, there is a long-range interaction between the fermions. 
Since the Lagrange multiplier field lives on three spatial and time dimensions, the interaction 
between fermions mediated by the exchange of the Lagrange multiplier fields is a logarithmic 
potential. Under the effect of such a confining potential, the fermions are confined. This is, of 
course, consistent with the constraint. Because the constraint is satisfied for low-energy exci- 
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tations. Contrary, if the fermions were massless, the propagation of the Lagrange multipher 
fields would excite many particle-hole pairs. The interaction between fermions mediated by the 
exchange of such a field is short range interaction. Therefore, there would be no confinement, 
or fermionic excitations would appear in the low-lying excitations, and the constraint would 
be no longer satisfied. This picture suggests that a small mass term in the Lagrange multiplier 
field should break the confinement of the fermions. This contradicts with the conclusion in 
ref.47 in which it is argued that the dynamical mass generation occurs in the presence of a 
small gauge field mass. 

As for the application to the high-temperature superconductors, the theory suggests that 
the low-lying excitations in the spin disordered regime is described by the fermions with the 
background of the 7r-fiux state correlations.^^ Because in the spin disordered regime there 
is no staggered magnetization. Therefore, there is no mass gap for particle-hole excitations. 
In this case the interaction between the fermions mediated by the Lagrange multiplier fields 
is of short range interaction. We may neglect the effect of such a short-range interaction. 
Meanwhile, the phase fluctuations are important for the properties of the fermions. Because 
the phase fluctuations play the role of transverse gauge fleld living in three spatial and time 
dimensions. As discussed in ref.49, such the transverse gauge fleld fluctuations lead to non- 
Fermi liquid behavior. 

In the continuum theory, or the QED3 theory, the gauge fleld associated with fluctuations 
about the mean fleld state is either compact or non-compact. By contrast, in our theory the 
phase fluctuations are compact but the Lagrange multiplier flelds are not compact. This may 
require a different approach to instanton effects discussed in the conflnement phenomenon in 
the compact QED3 theory.^'' 
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Appendix: Derivation of the Effective Interaction 

In deriving the effective interaction between the fermions, eq.(36), we first integrate over 
the fermion fields to obtain the effective action for the auxiliary fields, and then the effective 
interaction between the fermions is derived using the resulting effective action for the auxiliary 
fields. In this Appendix, we show that this kind of calculation is equivalent to the random 
phase approximation. To make clear the point, we consider a following simple model, 

H = JZ^kflfk + Y.<t>jUjk + Y,K,ct>,<t>-,. (A-1) 

k k,q q 
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Fig. A-1. The lowest order contribution to the fermion interaction due to the exchange of the (j>q 
field. Wavy line represents the propagator of the (j)q field and straight lines represent the fermion 
propagator. 






Fig. A-2. Feynnian diagram representing the random phase approximation. 



The lowest order term for the interaction between fermions mediated by the (pq field is given 

by 



VQ — 2 2^ J^ Jk+qJk'Jk'+qJk- 
k,k',q '^ 



(A-2) 



The Feynman diagram for this term is represented by Fig.A-1. By including fermion bubble 
diagrams and summing over terms represented by Fig. A-2, we obtain 



V, 



RPA 



lE 



1 



2 Z^ 2Ka + Ua(innr''+'^-"''-"'^'' 



Jk+aJk'Jk'+qJk, 



k,k' ,q 



i.q -r J-ig v""^'" 



where 



n, {i^n) = -^ E G'f i^^n) 4+, (i^n + i^n), 



I3N 



(A-3) 



(A-4) 



iUnM 



with GkiiuJn) = l/iiuJn - ^k)- 

Now we show that the same interaction term is obtained by the following path-integral 
calculation. We consider the action for the Hamiltonian (A-1), 



IUJ„ i—lUJ. 



Ik V-^^n 



n + Ck) fk {i(^n) +'^4>q (i^n) fl+q {i^n + i^n) fk {i^n) 
k,q 



+ E ^1^1 (^^") ^-<i (~^^n)- 



(A-5) 
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Integrating out fermion fields, we obtain 

S^ = -Tr In [{-iiOn + Ck) h,k'^n,n' + 4'k^k' {i^n - i^n')] + ^ KgCpg (iQn) 4>~q {-i^n) 



E 



Kg + -Uq {iUr 



(iQr. 



-i^n) + ■■■■ 



(A-6) 



Using the effective action for the (pq field obtained this way, we consider the following action, 



-iiOn + ^k) fk (i^n) + X] <Ag («^n) fl+q {i^n + i^n) fk (i^n) 
k,q 



+E 



Kq + -Uq (mn) 



<t>q (i^n) 4>-q {-i^n) ■ 



(A-7) 



By integrating out the (pq field, we obtain (A-3). 
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